Conductance of Atomic-Sized Lead Contacts in an Electrochemical Environment 
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Atomic-sized lead (Pb) contacts are deposited and dissolved in an electrochemical environment, 
and their transport properties are measured. Due to the electrochemical fabrication process, we 
obtain mechanically unstrained contacts and conductance histograms with sharply resolved, indi- 
vidual peaks. Charge transport calculations based on density functional theory (DFT) for various 
ideal Pb contact geometries are in good agreement with the experimental results. Depending on 
the atomic configuration, single-atom-wide contacts of one and the same metal yield very different 
conductance values. 

PACS numbers: 73.23.Ad, 73.63.Rt, 82.45.Yz, 31.15.es 



Charge transport through nanostructures presently 
constitutes a highly active area of research, also due to 
its relevance for the further miniaturization of electronic 
devices [T]. Regarding atomic-sized metallic contacts, 
different experimental techniques have been developed 
for their fabrication, e.g. mechanically controllable break 
junctions (MCBJs) [2_, modified scanning tunneling mi- 
croscopes (STMs) [3], electromigration [J], and electro- 
chemical methods [5] [3]. While in MCBJ experiments 
and STM-based setups atomic-sized contacts are formed 
by plastic deformation of the contact area, the advantage 
of the electrochemical deposition method is the absence 
of mechanical strain during contact formation, resulting 
in corresponding differences in conductance histograms. 

Conductance histograms of s-valent metals, such as 
gold, silver, or copper, have been analyzed extensively 
[TJ [7J showing a preference for multiples of Go = 
2e 2 /h, with the peak at IGq being a very robust fea- 
ture. As it is known that the chemical valence plays an 
important role for electron transport [TO] , conductance 
histograms for multivalent metals may consequently show 
a very different shape [T| |S] • 

In this work, Pb will serve as a model system for mul- 
tivalent metals. For Pb a broad peak between IGq and 
3Go was observed in conductance histograms based on 
the MCBJ technique [TT1 112]. During stretching of the 
junctions, the plateau regions of the conductance traces 
commonly exhibit a negative slope as a result of the 
mechanical strain US]. With electrochemical tech- 
niques we have recently fabricated atomic-scale silver 
contacts with a high thermal and mechanical stability 
[14] , and exploited this for constructing gate-controlled 
atomic switches [15] . In this work, we study charge trans- 
port through atomic-sized Pb contacts in experiment and 
theory. The conductance histogram is compared to the 
literature and theoretical results for ideal contact geome- 
tries. In the DFT-based calculations we explore both 
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Figure 1: (Color online) Schematic diagram of the experimen- 
tal setup. Pb is deposited and dissolved electrochemically in 
the narrow gap between two gold working electrodes (WEs) 
on a glass substrate, and the conductance of the contact is 
measured simultaneously. The potentials of the WEs with re- 
spect to the lead reference electrode (RE) and the lead counter 
electrode (CE) are set by a potentiostat. 



"single-atom" and "dimer" geometries with a single atom 
or a chain of two atoms in the narrowest part of the junc- 
tion, respectively. For these configurations we consider 
orientations of the semi-infinite electrodes along the three 
main crystallographic directions, namely (100), (110), 
and (111). 

Our experimental setup is illustrated schematically in 
Fig. [I] Since lead oxidizes under ambient conditions, the 
electrochemical cell is shielded in an inert gas chamber. 
Before starting with the electrochemical deposition, the 
chamber is thoroughly streamed with Ar. In this way, 
a practically oxygen-free environment is created for the 
growth of the atomic-sized Pb contacts. For their fabri- 
cation, two working electrodes (WEs) made of gold and 
separated by a gap of approximately 50 nm are prepared. 
They are covered with an insulating polymer coating ex- 
cept for the immediate contact area to minimize ionic 
conduction. Two lead wires (0.5 mm diameter, 99.998% 
purity) are used for the quasi-reference electrode (RE) 
and the counter electrode (CE). The potentials of the 
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Figure 2: (Color online) Conductance histogram of electro- 
chemically fabricated atomic-sized Pb contacts. The inset 
shows two typical conductance-time traces for deposition and 
dissolution processes. 



WEs with respect to the RE and the CE are set by a 
computer-controlled potentiostat. The electrolyte con- 
sists of ImM Pb(N0 3 ) 2 + 0.1 M HN0 3 in bi-distilled wa- 
ter. For conductance measurements an additional voltage 
of -12.9 mV is applied between the two WEs. To exclude 
possible specific effects of the gold WEs. contacts were 
also grown using Pb electrodes, yielding comparable re- 
sults. 

To fabricate the lead contact within the gap between 
the two WEs, a potential of 10-20 mV is applied to the 
RE. While lead is deposited in the junction, we monitor 
the conductance between the WEs. After contact depo- 
sition, the contact is dissolved again by setting the elec- 
trochemical potential of the RE to a value between -18 
mV and -36 mV, and closed again by choosing the poten- 
tial between 6 mV and 15 mV. By continuously repeating 
this procedure, conductance-time traces are recorded for 
a large number of opening and closing processes at a sam- 
pling rate of 50 ms. 

Two typical conductance-time traces obtained for con- 
tact closing (deposition) and opening (dissolution) are 
shown in the inset of Fig. [2j The plateau regions in the 
traces are quite flat in general. This fact can be explained 
by the lack of mechanical strain during the electrochem- 
ical growth process, in contrast to contacts fabricated 
with the MCBJ or STM techniques, where mechanical 
deformations are involved [HI EH H3J. The dissolution 
curve exhibits two plateaus at around 2.8Go and IAGq. 
These values are consistent with the range of the rather 
continuously decreasing last conductance plateau, as re- 
ported in Refs. [5] E2 The main panel of Fig. |2]shows the 
conductance histogram, plotted with a bin-size of 0.05Go- 
It represents the data from 1.5 • 10 6 conductance terraces 
in the range between and 20Gn for both deposition 
and dissolution from 64 different samples, each terrace 
being stable longer than 200 ms. The structure of the 
histogram with its most dominant peak at 1.4Go differs 
clearly from those obtained with the MCBJ approach, 



where only a single, broad peak was observed between 
IGo and 3Gn and no detailed structure within this peak 
could be resolved H~2] . Here, most probably due to the 
lack of mechanical strain and structural defects within 
the junction detailed substructure is resolved 

within this peak, demonstrating the advantages of the 
electrochemical fabrication method. The peak position 
at 1.4Go, however, is consistent with the results reported 
in Refs. ED E3 Additional peaks in the conductance dis- 
tribution are found at IGo, 2Go, and 2.8Go- In the range 
from 3Go to 6Go, the distribution is quite flat and shows 
a low broad maximum centered at around 4.8Gn. 

In order to understand better the experimental find- 
ings, we have performed a theoretical analysis of charge 
transport through various Pb contacts. We describe their 
electronic structure at the level of DFT and determine 
the conduction properties within the Landauer approach 
using Green's function techniques. Employing the RI- 
DFT module of TURBOMOLE 5.9 [TB], we use BP86 as 
exchange-correlation functional |17l [T5] and the standard 
basis set "def-SVP" of split-valence quality augmented 
with d-like polarization functions |TS] . An effective core 
potential efficiently deals with the innermost 78 electrons 
[2"U] . Our approach is described in detail in Ref. |2T| Let 
us note that we determine parameters for the description 
of the electrodes from a spherical fee cluster composed 
of 429 Pb atoms. The lattice constant for this cluster is 
set to the experimental value of 0.495 nm. The calcula- 
tion yields Ep = —3.76 eV for the Fermi energy, a value 
roughly corresponding to the experimental work function 
of 4.25 eV [22]. 

We investigate two types of junction configurations, 
namely single-atom and dimer contacts (Figs. [3]and|4|. 
For these we study different crystallographic orientations 
of the electrodes, namely (100), (110), and (111). The 
electrodes are oriented along the z axis, the direction of 
charge flow. Such geometries are believed to be responsi- 
ble for the first peak in conductance histograms of metals 
under ultra-high vacuum conditions. We consider only 
ideal geometries, where all atoms are located in the posi- 
tions of the fee lattice. The extended central clusters of 
Ref. I2T1 used for the description of the contacts, are com- 
posed of around 300 atoms. To keep the numerical effort 
manageable, their point group symmetry is exploited. 

Results for the single-atom contacts are displayed in 
Fig. [3] From the transmission at the Fermi energy, we 
obtain conductances of 2.7Go, 5.3Go, and 2.8Go for the 
(100), (110), and (111) directions, respectively. At first 
glance, the conductance for the (110) direction is sur- 
prisingly high. However, as is visible in Fig. |3jd), this 
contact should better be considered as a "five-atom" con- 
tact due to the additional bonds resulting from the small 
distance between atomic layers. Therefore, we will hence- 
forth refer to single-atom contacts as those for the (100) 
and (111) directions. For the latter two structures, the 
transmission at Ep is dominated by three transmission 
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Figure 3: (Color online) Geometries of single-atom contacts, their transmission as a function of the energy, and the local density 
of states (LDOS) of that atom in the narrowest part of the junction, which is indicated by an arrow in the geometries. The 
transmission r = is resolved into the contributions Tj from individual transmission channels, and the LDOS into its 

orbital components [21] . The crystallographic orientation of the semi-infinite electrodes is (a-c) (100), (d-f) (110), and (g-i) 
(111), respectively. Vertical dashed lines indicate the Fermi energy at Ef = —3.76 eV. 



channels. In each case, the non-degenerate channel is of 
sp z character, while the two degenerate ones are p x - and 
Py-like. The local density of states (LDOS) of the atom in 
the narrowest part of the constriction illustrates further 
the dominant role of the 6s and 6p states for conduction. 

Due to the directional character of transport the point 
group symmetries relevant for the classification of trans- 
mission channels in terms of irreducible representations 
are C^ v , C 2v , and G 3 „ for the structures with (100), (110), 
and (111) orientations, respectively, while their geomet- 
rical point groups are D4/1, D2h, and D^- Degeneracies 
of channels, present for the (100) and (111) configura- 
tions, are therefore removed for (110). This fact will be 
more clearly visible for the dimer junctions with their 
lower conductance, but the same point-group symme- 
tries, which we shall discuss next. 

In Fig. [4] we illustrate the evolution of the conductance 
and of the total energy as a function of the distance d be- 
tween the tip atoms for all the dimer contacts. All other 
distances are kept fixed in the stretching process. At the 
point of minimum total energy, conductances are 2.3Go, 
I.6G0, and 2.8Go for the (100), (110), and (111) con- 
tacts, respectively. As naively expected, the conductance 
exhibits a monotonous decay with d for (100) and (110). 
For (111), however, it stays rather constant. From an 
analysis of the LDOS of one of the chain atoms, we find 
that this is due to the p x and p y states, which move into 
resonance. As mentioned before, degeneracies of con- 
duction channels are lifted for the (110) orientation as 
compared to the other junction configurations. 



The peaks at 1.4Go and 2.8Go in our histogram of 
Fig. [2] are consistent with our theoretical results for 
single-atom and dimer structures. This indicates that 
the contacts observed in the experiment resemble ideal 
lead junctions without strong contaminations or struc- 
tural defects, since they would lead to deviations from 
these conductance values. 

The experimentally observed, pronounced peak at IGo 
is not predicted by the conductance calculations de- 
scribed above. However, as the contacts are in an aque- 
ous electrochemical environment, small molecular species 
working as bridges between the two atomic-scale Pb elec- 
trodes, such as traces of H2, O, O2, or H2O, might play 
a role |12l |2"51 |2"4"] . Alternatively, ionic lead species could 
be taken into account, where specific transport channels 
of the junction are selectively suppressed. The interpre- 
tation of the peak at IGo as a result of the conductance 
quantization of a free electron gas in a geometrical con- 
striction, as observed in semiconductor heterostructures 
[251 126] and s-type metals pQ, appears as an oversimpli- 
fication for multivalent metals |12[ 127], 

An alternative explanation of the peaks at 2Go and 
2.8Go would be their interpretation as replicas of those 
at IGo and 1.4Go with twice the value. Although geome- 
tries with parallel atomic wires were observed for STM- 
fabricated gold contacts [21]) no indications for the for- 
mation of similar structures exist for lead so far. 

Our analysis of transport using DFT is in agreement 
with theoretical results reported in Ref. 1131 Employing a 
simpler tight-binding parametrization for describing the 
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Figure 4: (Color online) Dependence of the conductance G and the total energy Etot on the distance d between the tip atoms 
for the different dimer configurations with (a-b) (100), (c-d) (110), and (e-f) (Iff) crystallographic orientation. In addition to 
G = Gi the contributions from the individual transmission channels Gi = Gon are shown. In each plot the y-axis on the 
left refers to G and Gi (solid lines) and that on the right to E to t (dashed lines). Dotted vertical lines indicate the distances d 
with minimal E to t- 



electronic structure and considering only the single-atom 
contact oriented along the (111) direction, the authors 
report the same conductance of 2.8Go- Keeping the cen- 
tral atom fixed, they find that the conductance stays 
practically constant when the junction is stretched, sim- 
ilar to our observations for the (111) dimer structure. 
They argue that the negative slopes of the conductance 
plateaus before contact rupture would require the inclu- 
sion of spin-orbit interactions. While such relativistic 
effects are considered in our calculations only in an effec- 
tive manner through the core potentials [20], our results 
show that the dependence of the conductance on distance 
is strongly geometry-dependent. 

To conclude, the high thermal and mechanical stability 
and low strain in our clcctrochcmically fabricated atomic- 
sized lead contacts makes it possible to resolve sharp 
conductance peaks and detailed substructures within the 
broad first peak in the conductance histograms obtained 
with MCBJ or STM techniques. The theoretical analysis 
of transport within our DFT approach yields junction ge- 
ometries compatible with the peaks at 1.4G*o and 2.8G*o 
and allows to attribute them to neutral single-atom and 
dimer contacts of Pb. The results demonstrate that the 
conductance of Pb contacts with a single-atom- wide con- 
striction depends crucially on the contact geometry. 
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